Selective Mottness as a key to iron superconductors 
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The discovery of high-T c superconductivity in 
iron pnictides[ ] has opened a new path to- 
wards the understanding of this phenomenon. In 
the cuprates, the materials with the highest T c , 
electron-electron interactions are the unifying el- 
ement between the anomalies of the metallic state 
to the mechanism of superconductivity. In the 
iron family a likely electronic (magnetic) super- 
conducting pairing [ ] is not accompanied by a dis- 
tinctive evidence of strong correlations. 

Here we show that the experimental phase di- 
agram of electron- and hole-doped BaFe2As2 can 
be unified - the stochiometric valence of 6 elec- 
trons per iron playing no special role for corre- 
lations - and that these materials are influenced 
by a Mott insulating state that would be real- 
ized at a doping of one hole/Fe (i.e. for half-filled 
conduction bands). The key factor is a strong 
Hund's coupling, which decouples the different 
orbitals. Each orbital behaves as a single band 
model, where the correlation degree only depends 
on how doped is each orbital from half-filling. As 
a result weakly and strongly correlated conduc- 
tion electrons coexist in much of the phase dia- 
gram, an effect which increases with hole doping. 
This reconciles contrasting evidences on the elec- 
tronic correlation strength and places the super- 
conducting phase around 20% average doping per 
orbital from a half-filled Mott insulator, as in the 
cuprates. 

Electronic correlation induces in solids a wide palette 
of unconventional behaviors, not captured by the stan- 
dard Density Functional Theory (DFT), where electron- 
electron interactions are treated at the one-body level. 
A spectacular example is high-Tc superconductivity in 
the cuprates, which is broadly believed to origin from the 
same correlations that make the parent compounds Mott 
insulators. In this light the question of the correlation 
strength in iron superconductors, the runner-up materi- 
als in terms of T c (the highest to date being 57K), is one 
to be clarified. Contrasting experimental evidences have 
been reported, and early theoretical approaches start- 
ing both from the weak-coupling [ ] and from the strong- 
coupling[ ] limit have been implemented. 

Indeed despite a bad metallic conductivity and many 
other unconventional properties [4, 5], experiments have 
failed to detect a direct signature of strong local inter- 



actions such as a Hubbard band[ 6] and have instead re- 
ported an overall fair agreement with DFT predictions 
and assessed the local screened Coulomb repulsion U in 
the d-orbitals to be smaller than the total width of the 
conduction bands W ~ 4eV. 




Figure 1. Experimental mass enhancements in doped 
BaFe2As2 assessed by different techniques (see online supple- 
mentary information for references). The experimental phase 
diagram, including the antiferromagnetically ordered metallic 
phase (AF) and the superconducting one (SC) is reproduced 
as a background. The solid line is a fit of the specific heat 
data. ARPES and quantum oscillation data corresponding to 
a given doping represent the values estimated for the differ- 
ent sheets of the Fermi surface. With reducing filling correla- 
tions increase monotonously, and the estimates from different 
techniques spread more and more, an indication of increasing 
selective Mottness (see text). 

However it was emphasized early on[7] that Hund's 
coupling J plays a major role in these compounds. 
This term, which favors high-spin configurations on each 
atom, is consistently estimated to be ~ 0.7 — 0.9eV both 
theoretically and experimentally. On theoretical grounds 
it is shown[<s] that for the filling of 6 electrons in 5 orbitals 
it favors a bad-metallic behavior yet pushing away the 
Mott insulator. It has also been shown [7] that this kind 
of correlated metal does not display as prominent Hub- 
bard bands as in standard correlated materials. Overall 
an image of moderately strong correlations [9] driven by 
Hund's coupling, seems to emerge. 

In the present paper we take a slightly different point 
of view, and show how in most of the phase dia- 
gram of the T22-family' (i.e. both hole- and electron- 
doped BaFe2As2) the degree of correlation is heavily 
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differentiated among the electrons forming the conduc- 
tions bands. This brings about a scenario of coexist- 
ing strongly and weakly correlated electrons. Indeed 
it was shown by two of us that iron superconductors 
are the ideal ground for heavily orbital-dependent cor- 
relation effects (selective Mottness), or even for orbital- 
selective Mott localization[l(J, 11]. Phenomenological de- 
scriptions assuming the presence of itinerant and local- 
ized electrons [12- 14] have indeed reproduced many fea- 
tures of these materials. 

In Fig. 1 we have collected experimental mass enhance- 
ment estimates from specific heat, infrared reflectivity, 
angle-resolved photoemission spectroscopy (ARPES) and 
quantum oscillations (the references are detailed in the 
supplementary online information, along with a discus- 
sion of the reliability of the estimates) . All the data are 
taken in the high-temperature tetragonal metallic phase 
above the magnetically ordered (and orthorombically dis- 
torted) phase at low dopings and the superconducting 
one at intermediate dopings. 

Single-crystal specific heat measurements delineate a 
very clear trend: the mass enhancement m* /nib = 7/76 
(where 7 and 75 are respectively the measured normal- 
state Sommerfeld coefficient and the band theory esti- 
mate 7b = ir 2 kg/3No, No being the total density of 
states at the Fermi energy) grows continuously from ~ 2 
in the electron-overdoped Ba(Fei_ :r Co :E )2As2 to 3 in the 
electron-underdoped/stochiometric compounds towards 
5-6 in the optimally hole-doped Bai_ x K a: Fe2As2 until 
reaching the very high value of ~ 9 for the fully sub- 
stituted KFe2As2 (having a doping of 0.5 holes/Fe). 

Optical conductivity provides a partially different pic- 
ture. Comparing the band-theory estimated of the Drude 
weight 

D band = J2(g^) 2 m-^), (1) 

(where ej? is the dispersion of band a) with the ex- 
perimental value one obtains an optical effective mass 
Dband/ D exp = m* /mi,. With this approach Qazilbash et 
al.[13] found intermediate correlation strength (~ 3.3 for 
BaFe2As2 while the data we collect show correlations di- 
minishing with electron doping but a possible saturation 
(with m* I rrib hardly exceeding ~ 4) of the mass enhance- 
ment for the hole-doped compounds. 

A general argument links this apparent discrepancy to 
the coexistence of electrons with different effective mass 
renormalization: in this multiband environment both the 
conductivity and the density of states at the Fermi en- 
ergy stem from a sum over the band (or orbital) index. 
As a consequence 7/76 = Db an d/ D exp only if the mass 
enhancement is the same for all electrons. If instead the 
renormalization is different from band to band (from or- 
bital to orbital) the density of states is a sum of contri- 
butions ~ (m* /nib) a whereas the Drude weight is a sum 



of terms ~ (m&/m*) a , and thus as it happens for series 
or parallels of resistances, the first will be dominated by 
the largest mass enhancement, the second by the small- 
est. This means that in case of heavy differentiation of 
correlation strength the enhancement of the Sommerfeld 
coefficient will reflect that of the strongly correlated elec- 
trons, whereas the Drude weight the one of the weakly 
correlated ones. 

A more direct confirmation of this picture is obtained 
by means of ARPES. In order to roughly match the ex- 
perimental data the DFT bandstructure has to be renor- 
malized by a global factor (indicated as 'ARPES whole' 
in Fig. 1) that increases with decreasing filling. More- 
over by measuring the Fermi velocity on the different 
Fermi sheets and comparing them to DFT calculations 
one can estimate independently the mass enhancements 
of the corresponding low energy electrons. As we report 
in Fig.l, in the electron-doped side of the phase diagram 
the different Fermi sheets show a homogeneous renor- 
malization of the Fermi velocities by a factor of order 
~ 2, whereas moving to the hole-doped side the different 
Fermi sheets show different renormalization for different 
bands (between 3 and 9 for Bao.6Ko.4Fe2As2 and between 
2 and 18 for KFe2As2). Quantum oscillations measure- 
ments performed on KFe2As2 also report a heavy differ- 
entiation between the mass enhancement. 

Thus in summary the unified phase diagram of 
electron- and hole-doped BaFe2As2 shows a strong and 
monotonous increase of correlations going from the heavy 
electron-overdoping to the heavy hole-overdoping, and a 
parallel progressive differentiation of the mass enhance- 
ments, or selective Mottness. It is also worth noting that 
the stochiometric filling does not play a particular role 
in this phase diagram of the tetragonal phase, contrary 
to the fact that commensuration is naively expected to 
strengthen the effects of correlations. 

Calculations including dynamical correlations in DFT 
ab-initio calculations for both iron pnictides and chalco- 
genides show that the combined role of U and J gives 
rise to a novel scenario (as summarized in Ref. [16]). A 
sharp crossover[17-19] is found, as a function of increas- 
ing correlation strength and decreasing filling, towards a 
'Hund's metal' [20] state. The physics there is dominated 
by Hund's coupling and can be described using a few con- 
cepts isolated in model studies: 1) a strongly correlated 
metallic phase far from the Mott transition at n=6 [8] is 
realized[18, 19, 21] with a very low coherence tempera- 
ture and correlation strength more sensitive to J than to 
U[7]; 2) above such coherence temperature a 'spin-frozen' 
state appears[22, 23], in which the self-energy has a pe- 
culiar y/uj dependence at low u and long spin correlation 
time[24]; 3) 'orbital decoupling' [ ] happens i.e. differ- 
entiation of the correlation strength among the different 
orbitals, due to the fact the Hund's coupling suppresses 
interorbital correlations, rendering charge excitations in 
the various orbitals virtually independent from one an- 
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Figure 2. Calculated orbitally resolved mass enhancement 
(m*/m(,) a and quasiparticle weight Z a (two uppermost 
panels) and electron population (middle panel) for doped 
BaFe2As2 and stochiometric KFe2As2 (red dots) within Slave- 
Spins mean-field. Note that for each orbital a, Z a tends to 
vanish when n a approaches half-filling. Indeed when plot- 
ting Z a as a function of n a for each orbital (bottom panel) 
a striking linear behavior appears showing that each orbital 
is only sensitive to its doping away from the half-filled Mott 
insulator found at total population n=5, independently from 
the others. The dots signal the values for the stochiometric 
compound, having total population n=6. 



other. 

This last point is what we highlight here as the main 
reason for the experimental scenario we have detailed. 

We have performed LDA+Slave-spin mean-field (see 
'Methods' section) calculations at T = starting from 
DFT bandstructures for doped BaFe2As2 in the tetrago- 
nal phase (see the Supplementary Information for other 
iron superconductors). We have checked that for the 
interaction strength corresponding to the ab-initio esti- 
mate this compound displays the signature of the 'Hund's 
metal' phase. Our results as a function of U (see supple- 
mentary information) show that BaFe2As2 lies around or 
just above a cross-over to a correlated phase signaled by 
a fast drop in the quasiparticle weight Z, an increase of 
the interorbital spin correlation functions and decrease in 
the inter-orbital charge correlation functions, as reported 
in [18, 19]. As expected in this regime the correlation 
strengths are strongly differentiated among orbitals, with 
the tig orbitals more correlated the e g ones; in particular 
the xy orbital is the most correlated of all. This mo- 
tivates the experimental observation from ARPES that 
at low energy (where bands are mainly of ti g character) 
the bandstructure is more renormalized than at high en- 
ergy (where both i 2 g and e g characters are present). This 
heavy-fcrmionic behavior might be turned into an inco- 
herent phase at T =/= and/or with a refined treatment 
of quantum fluctuations[22, 23]. 

We have studied the evolution as a function of to- 
tal density of the orbitally resolved quasiparticle weights 
Z a = (m\ 1 /m*) a and populations, which are shown in 
Fig. 2. Correlations clearly increase continuously de- 
creasing the total population, resulting in a more cor- 
related hole-doped side compared to the electron-doped 
one[2G]. The commensurate filling of n=6 does not play 
any special role in this undistorted phase. 

It can be noted that the theoretical mass enhancements 
reach values, for strong hole overdoping, much higher 
than experiments. This is explained by a progressive in- 
crease of covalency in the Fe-As planes introduced by the 
Ba — » K substitution: when this is taken into account 
(using the DFT bandstructure for KFe 2 As2 at n=5.5) 
correlation strengths close to the experimental values are 
obtained (red dots in Fig. 2). 

When the suggested 'orbital decoupling' mechanism is 
realized the system switches from a collective to an indi- 
vidual orbital behavior: the correlation strength in each 
orbital depends mainly on the doping of the orbital away 
from half-filling. We have tested this hypothesis here and 
found that it is fully verified: in the lower panel of Fig. 2 
we plot the quasiparticle weight of each orbital as a func- 
tion of its population. Each orbital is only sensitive to 
its own filling and behaves like a doped single-band Mott 
insulator, with Z proportional to the doping. With de- 
creasing total population all orbitals move towards half- 
filling and the quasiparticle weights go to zero at n=5, 
where the system is then a Mott insulator, with vanish- 
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Figure 3. Possible unified phase diagram for pnictides 
and cuprates. We plot here the experimental phase dia- 
gram for doped BaFe2As2 as a function of the average or- 
bital doping. In the hypothesis that the magnetically or- 
dered/orthorombically distorted phase (grey area) is acciden- 
tally favored by the total commensurate filling n=6 around 
the stochiometric compound thus suppressing the supercon- 
ductivity of the tetragonal phase, we artificially eliminate it 
and complete the superconducting dome (blue area). We then 
see that this dome is centered around 20% doping (per orbital) 
away from the half-filled Mott insulator (here at n=5), as in 
cuprates. In between we find a strongly correlated Mott se- 
lective (and bad metallic [17]) phase, while at higher orbital 
dopings (corresponding respectively to the electron doped re- 
gion in pnictides and the overdoped region in cuprates) we 
recover a moderately correlated Fermi liquid. 



ing orbital polarization. The remarkable thing is that, 
albeit this happens at a very different rate for each or- 
bital (the xy orbital is almost half-filled already at n=6, 
and thus more strongly correlated, the z 1 is quite far from 
half-filling until very close to n=5, and thus weakly cor- 
related for most of the fillings) the behavior when plotted 
against the individual population is quite universal (the 
different slopes being determined by the details of the 
bandstructure) . 

This scenario confirms what observed in previous the- 
oretical studies[17, 27], that the parent compounds for 
pnictides at n=6, are influenced by correlations due to the 
Mott insulator at n=5. Here we provide the reason why 
: the orbitals are 'decoupled' from one another thus ren- 
dering the individual population of each orbital crucial 
for its correlation strength. Accordingly, the seemingly 
far in population n=6 compound, has orbitals instead 
very close to the Mott insulating state. The ti g orbitals 
and the xy in particular are only a few percents of doping 
away from the selective Mott insulator. The coexistence 
of electrons whose different correlation degree descends 
from a different distance to the Mott insulator defines 
the selective Mottness and its role in the phase diagram 
This scenario is also consistent with the gap (which is 
the remnant of the n=5 Mott Gap) at negative energies 
in the spectral functions in Hartree-Fock mean- field [ ]. 



One last more speculative hypothesis is now in order. 
By plotting the experimental phase diagram as a function 
of the average orbital population the stochiometric filling 
falls at 20% average doping [I 7]. This is now justified by 
the orbital-decoupled nature of the conduction electrons. 

If, as quite generally believed, the magnetically or- 
dered (and orthorombically distorted) phase is detrimen- 
tal to superconductivity [4] and it is favored accidentally 
by nesting and/or lack of disorder realized for the com- 
mensurate n=6 filling, we can imagine that the uni- 
fied phase diagram of the tetragonal phase of iron su- 
perconductors would see the two domes unified into a 
bigger one (see Fig. 3). This tentative phase diagram 
would now be strikingly similar to that of cuprate su- 
perconductors: a superconducting dome would be cen- 
tered around a doping of 20% per band away from a 
Mott insulator, and flanked at smaller doping by a bad 
metallic phase in which electrons with different correla- 
tion strength coexist, as it happens in hot /cold spots in 
the underdoped cuprates - notice that momentum-space 
selectivity observed in cluster DMFT studies of the two- 
dimensional Hubbard model can be easily recast into a 
multi-orbital effective picture reminding the situation in 
iron pnictides[29]. At higher fillings a weakly correlated 
good metal is recovered in both cases. This scenario uni- 
fies a series a seemingly contradictory evidences and it 
is expected to put serious constraints on pairing theories 
for iron-superconductors. 
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EXPERIMENTAL ESTIMATES OF THE 
CORRELATION 

We report here the values and references of the exper- 
imental data shown in Fig. 1 of the main text, where we 
plot several estimates of the quasiparticle mass enhance- 
ment in the tetragonal phase of the 122 pnictide family. 
We also outline a few caveats on the data analysis of the 
various experimental probes, and we specify the physical 
significance of the mass enhancement extracted with each 
method. The original references are cited in the tables, 
and most of these results are nicely summarized in the 
reviews Refs. [1, 2]. 

Specific Heat 

Let's first discuss the mass enhancement as estimated 
from the specific heat. At low temperatures the specific 
heat of a Fermi-liquid goes like 

C v =-/T + 0(T 3 ) (2) 

where the "Sommerfeld coefficient" 7 is proportional to 
the renormalized mass of the carriers. The renormaliza- 
tion factor is then usually estimated as m*/m = 7/7)1, , 
where -fb = K 2 k 2 B /?>N{e F ) and N(ep) is the total non- 
interacting band density of states at the Fermi energy 
cf, as estimated by DFT. 

The experimental values for 7 used to calculate the 
mass enhancements plotted in Fig.l of the main text and 
the corresponding DFT values, along with the original 
references are reported in table SI. 

A source of error on 7 is however that in most of the 
cases reported here the specific heat of the normal tetrag- 
onal phase has to be extrapolated to zero temperature. 
Indeed the actual ground state is either a superconductor, 
or an orthogonally distorted phase with magnetic order, 
that we are not treating here. What is done then in or- 
der to estimate 7 of the high-T tetragonal phase is to fit 
Cy/T in the high-temperature phase and extrapolate it 
to zero T. When the low-T phase is superconducting this 
fit can be helped by the additional constraint due to the 
second-order nature of the superconducting transition, 
that imposes matching entropies S(T) — J dTCy /T at 
T c . Then from the actual measurement of Cy in the su- 
perconducting phase one can impose a constraint on the 
integral of Cy /T of the normal phase in the extrapolated 
range of temperatures. No similar constraint can be used 
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Table SI. Mass enhancements estimated from specific heat 
measurements used in Fig. 1 of the main text. The columns 
report respectively the values for the Sommerfeld coefficient 
7 extracted from experiments, the corresponding DFT calcu- 
lated values 7b and the calculated mass enhancement. Several 
LDA values for 71, are available for KFe2As2 in the literature 
ranging from 10.1 [8, 9] to 13.0[9]. Our calculation yields 14.5, 
thus confirming the wide error bar associated to this value. 
The estimated mass-enhancement ranges then form 6.5 to 9.3. 
We use the more dramatic value of 9.3 in the Fig. 1 of the 
main article in order to better illustrate this effect, but all 
values are higher than the ones estimated anywhere else in 
the phase diagram, thus confirming the solid physical claim 
of increasing correlations estimated with specific heat, with 
reducing total population. 

around the stochiometric n = 6 compound, where the 
system undergoes a first-order magnetic/structural tran- 
sition at > 100K, thus making the extrapolation quite 
unreliable. For this reason we omitted the data in this 
region. It is worth reminding that the reliability of these 
estimates are also subject to the ability to subtract the 
phononic contribution to the specific heat, that has to be 
estimated independently 

Also, the purity of the sample is of utter impor- 
tance. Only measurements performed on single crys- 
tals are shown in Fig.l of the main text, specific heat 
in polycrystals being subject to substantial variations. 
An example is the substantially reduced value 7 = 
69mJ/molK 2 reported[12] for polycrystalline KFe2As2 
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Figure SI. Sommerfeld coefficient 7 for single-crystals used in 
Fig.l of the main text, compared with the values measured in 
polycrystals. The gray area represents the zone where the low- 
T phase is orthorombical distorted and magnetically ordered, 
and where the specific heat of the high-T tetragonal param- 
agnetic phase has to be extrapolated from above the SDW 
transition. Outside this zone, where the error bars on 7 are 
necessarily higher, even polycrystals show an enhancement of 
the Sommerfeld coefficient with hole doping. Polycrystal data 
are taken from [10-12] 

of residual resistivity ratio RRR=67, to be compared 
with 7 = 94mJ/molK 2 reported [7] for a single crystal 
with RRR=650. See Ref. [1] for a somewhat more thor- 
ough discussion on this point. This being said, where 
their estimate is more reliable also polycrystals show an 
increasing Sommerfeld coefficient with hole doping (see 
Fig. SI). 

Optical conductivity 

Other mass enhancements reported in Fig.l of the 
main text are extracted from optical conductivity. In- 
deed the peak in the spectrum at low frequency (the 
Drude contribution) is usually considered to account for 
the quasiparticle contribution to the conductivity. Its 
width accounts for the scattering rate of these charge 
carriers, whereas its total weight for their effective mass. 
Comparing this weight with the value obtained in DFT 
is a possible way to assess the mass enhancement due to 
dynamical correlations. Indeed correlations reduce the 
Drude weight, and in a Mott insulator its value vanishes. 

This can be seen as an analysis of the reduction of 
the kinetic energy of the quasiparticles[13] due to cor- 
relations, K exp / 'Kb an d, that indeed vanishes in a Mott 
insulator. The overall kinetic energy of the conduction 
electrons however, is reduced by correlations but does not 
vanish, even in a Mott insulator where itinerant (albeit 
incoherent) excited states still exist at higher energy. The 
spectral weight of intraband transitions always adds up 
to the conduction electrons kinetic energy[ 14] but in cor- 
related materials it is transferred away from the Drude 
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Table S2. Mass enhancement values extracted from optical 
conductivity measurements and used in Fig. 1 of the main 
text. The mass enhancement is estimated via the ratio of the 
Drude spectral weight and the Drude+MIR spectral weight. 
Two method are used to estimate these spectral weigths, as 
described in the text. The first (results in columns 2) uses 
a fit to Drude-Lorentz oscillators, the second (column 3, the 
hole-doped values are integrals calculated by digitizing the 
data reported in the cited references) uses two cutoffs: £l c = 
500cm -1 for the Drude component, and fi c = 2000cm -1 for 
the Drude+MIR weight. The values of column 3 are used in 
Fig.l of the main text. 

peak into mid-infrared (MIR) and higher energy features. 
Thus an alternative, approximate way to estimate the rel- 
ative reduction of the quasiparticle kinetic energy based 
exclusively on experimental data was proposed by Lu- 
carelli et al.[15, 16]. It consists in normalizing the spec- 
tral weight of the Drude peak to the Drude+MIR total 
spectral weight. The latter is indeed an underestimate of 
the total band-theory kinetic energy, but has the advan- 
tage of excluding inter-band transition that contaminate 
the spectral weight at still higher energies. 

In Ref. [16] it was shown that mass enhancements esti- 
mated through this procedure agree rather well with the 
estimate K exp / K^and in BaFe 2 As2 of Ref. [13]. Thus this 
analysis was performed on the electron-doped compounds 
in [15] using two methods: the Drude and MIR contribu- 
tion were determined using a Drude-Lorentz model and 
then integrated separately, or they were approximately 
separated using two appropriately chosen cutoffs. The 
mass enhancement was defined as: 

TO* /m ~ Kband/ K exp ~ SWorude+MIB./ SW Drude, (3) 

with SW{Vt c ) = J Q c duitj'{u), where <r'(w) is the real part 
of the optical conductivity, fl c — 500cto -1 for the Drude 
contribution and £l c — 2000cto -1 for Drude + MIR. 

Here we extend this last method to available experi- 
mental data on the hole-doped side and report the ob- 
tained mass enhancements in table S2. 

It is worth pointing out that several technical points 
make this analysis highly nontrivial and are possible 
sources of error. It is difficult for instance to isolate 
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the Drude contribution within the optical conductivity: 
indeed interband transitions partially superpose to it, 
especially when the best structure is as intricate as in 
iron superconductors, where multiple low-energy inter- 
band transitions are available. This makes the choice of 
a cutoff for the spectral weight integral uneasy. On the 
other hand fits based on the Drude-Lorentz model can 
be used to sort out the Drude contribution, but often 
several alternatives are possible. An ubiquitous compli- 
cation is the need of at least two Drude components in 
order to fit the low-energy spectral weight. One Drude 
peak is narrow while the other is very wide, accounting 
for incoherent processes and overlapping with interband 
contributions. Last but not least, the coupling to bosonic 
modes (as magnetic fluctuations for example) can also 
shift some spectral weight, rendering the fits or the choice 
of the cutoffs even harder. 

Overall then, the mass enhancements estimated 
through the analysis of the optical spectral weight that 
we have performed have to be considered as indicative. 
And indeed they yield a picture of intermediate corre- 
lation strengths throughout the phase diagram of doped 
BaFe 2 As2, turning to rather weaker correlations in the 
electron-overdoped region [20]. 



ARPES and Quantum Oscillations 
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Table S3. Mass enhancement values extracted from ARPES 
and Quantum Oscillation measurements and used in Fig. 1 of 
the main text. What we highlight in this paper is that they 
spread moving from electron-doping to hole-doping. The first 
column reports the rough proportionality factor in order for 
the measured bandstructure and the DFT one to best match 
each other globally (showing that also the degree of correla- 
tion increases going from electron- to hole-doping. The follow- 
ing columns report the mass enhancements estimated as the 
ratio between the theoretical and measured Fermi velocities, 
resolved for every Fermi sheet. a,/3 and £ are hole-like pock- 
ets, whereas 7, S and e are electron-like ones. The electron 
doped data are taken at k z ~ 1 ((i.e. along Z — A) whereas the 
hole-doped are taken at k z ~ (i.e. along V — M. The last 
column reports the mass enhancement of the band of main 
character 3z 2 — r 2 that lies below the Fermi level and slightly 
crosses it near Z in KFe2As2. Finally Quantum Oscillations 
measurements (that do not observe the ft Fermi sheet) for 
KFe2As2 are reported. 



We report in table S3 the mass enhancements ex- 
tracted from ARPES measurements in both electron- and 
hole-doped BaFe2As2 and in KFe2As2, for which we also 
report Quantum Oscillation measurements. 

What we highlight is that going from electron-doping 
to hole-doping until reaching KFe2As2 (0.5 holes/Fe dop- 
ing) correlation strength increases (looking at the global 
renormalization factor needed in order to match the theo- 
retical and the observed bandstructure) . These data also 
show the measured mass enhancement for each Fermi 
sheet (each of which bears a different dominant orbital 
character), spread more and more with the reduction of 
the filling, thus confirming the increasing orbital selec- 
tivity. 

A caveat is due however when directly comparing 
ARPES and quantum oscillations with the DFT band- 
structures. Electronic correlations induce both band- 
width renormalizations and band shifts. This can result 
in a change of the location (and size) of the Fermi surface 
and thus comparing the electron velocities at the Fermi 
level can be misleading. In principle a safer procedure 
would be to compare the measured Fermi velocity with 
the calculated velocity of bands once shifted in order to 
better match experiments] ]. Despite the sizable error 
bars induced by this effect, the large spread of the mass 
renormalizations depending on the different character of 
electrons at the Fermi level is clear trend. 



Other experimental data supporting selective 
Mottness in Iron Superconductors 

Besides some of the cited studies (the optics papers 
(Lucarelli et al.[15] and Wang et al.[19]) interpreting re- 
sults as a coexistence of more itinerant and more localized 
electrons and the ARPES studies by Ding et al.[23] and 
Yoshida et al.[24] highlighting the strong orbital depen- 
dence of correlations) other experiments point towards a 
main role played by selective Mottness in Iron pnictides. 
Yuan et al. [20] argue that magnetotransport data in 
Co-doped BaFe2As2 can be interpreted as showing that 
both itinerant and localized electrons are present, and 
that the scattering of the former on the latter gives rise 
to the complex transport properties that are measured. 
Malaeb et al.[27] have shown a strong orbital-dependence 
of the superconducting gap with hole-doping, leading 
to a disappearance of only one gap in superconductive 
Ba .4Ko.6Fe 2 As2 

Also in the arguably more correlated Iron chalco- 
genides several experiments pointed towards a localized- 
itinerant dichotomy. NMR and EPR measurements on 
FeSeo.42Teo.58 by Arcon et al. [28], and neutron scat- 
tering measurement on FeSeo.65Teo.35 by Xu at al. [ | 
are intepreted as showing the presence of intrinsic local 
magnetic moments in the metallic non-superconducting 
phase, coexisting then with itinerant electrons. Also Yi et 
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al. [ )] study the intercalated chalcogenides A a; Fe2-j / Se2 
(A=K, Rb) and report the disappearance of the band 



of xy-character when the temperature is increased above 
~ 150-ftT, thus signaling the orbital-selective Mott transi- 
tion of the most correlated electrons. 
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Figure S2. Orbitally resolved quasiparticle weights and Spin-spin and charge-charge interorbital correlation functions for 
BaFe2As2, as a function of U for J/U=0.25. 
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Figure S4. Orbitally resolved quasiparticle weights and Spin-spin and charge-charge interorbital correlation functions for FeSe, 
as a function of U for J/U=0.224. 



THEORY 



Methods 



The material-specihc band structures are determined 
using Density Functional Theory with the Generalized 
Gradient Approximation for the exchange-correlation po- 
tential according to the Perdew-Burke-Ernzerhof recipe 
as implemented in Quantum Espresso [31]. Then we 
apply Wannier90 [32] to compute the maximally local- 
ized Wannier orbitals, which we use as a basis to gen- 
erate a tight-binding Hamiltonian for the 5 orbitals of 
d-symmetry giving rise to the conduction bands. We in- 



clude the electron-electron on-site interaction in the form 
H =U ^ n imt n imi + U' ^ n ima n lmS (4) 

i,m i,m>m f ,cr 
"F (U </) ^ ^ ^iraa^im' (j 

where Uj mCT = d| m<T di mff and di ma is the destruction op- 
erator of an electron of spin a at site i in orbital m. U 
and U' = U — 2 J are intra- and inter-orbital repulsions 
and J is the Hund's coupling, treated here in the Ising 
(density-density only) form. 

This 5-orbital Hubbard-like Hamiltonian is solved 
within the Slave-Spin mean-field scheme, here general- 
ized to the case of realistic hamiltonians. The original 
implementation of this method[33, 34] (and differently 
from a more recent implementation [ 35]) yields nonzero 
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Lagrange multipliers in the noninteracting limit, when a 
crystal-field is present. A constant was added to each 
Lagrange multiplier in order to compensate this artifi- 
cial level shift, and adjusted in order to obtain the cor- 
rect non-interacting population distribution among the 
orbitals at zero interaction strength. 

Details of the implementation will be given 
elsewhere I'M)]. 



compound. Indeed, as illustrated in Figs S5 and S6 both 
LaFeAsO and FeSe show a very similar behaviour to the 
one reported for BaFe2As2 in Fig. 2 of the main article. 
The orbitally resolved quasiparticle weights are linearly 
dependent on the population of each orbital. Indeed cor- 
relations increase when reducing the total population, the 
system approaching the n=5 Mott insulator. 



Orbital-decoupling in other iron superconductors 

In this section we compare slave-spin results for sto- 
chiometric and doped BaFe2As2, LaFeAsO and FeSe. 

Given the approximate nature of our many-body treat- 
ment (for instance a somewhat larger J is needed in 
the present formulation of the slave-spin mean-field com- 
pared to DMFT) , we did not use the exact ab-initio es- 
timates of the interaction parameters. For the Hubbard 
repulsion we use U = 2.7eV for BaFc 2 As 2 , U — 2.7eV 
for LaFeAsO U = A.leV for FeSe, which are also close 
to cRPA estimates for a d-model[37], whereas we explore 
a range of values for the Hund's coupling J. For the cal- 
culations on BaFe2As2 shown in Figure 2 of the main 
article we have chosen J/U = 0.25 in order to reasonably 
reproduce the experimental effective mass at n=6. We 
then use J/U = 0.29 for LaFeAsO and J/U = 0.224 for 
FeSe. 

We first present the results for the stochiometric com- 
pounds as a function of U, keeping the ratio J/U fixed. 

These data show a striking similarity between the dif- 
ferent compounds: after an initial slow reduction of the 
quasiparticle weights Z at weak interaction strength, they 
all undergo a crossover in which the quasiparticle weights 
quickly drop and the system enters a different metallic 
phase (the Hund's metal). Here the Z are essentially 
flat and different among the orbitals. The spin-spin cor- 
relation functions saturate to high values, indicating the 
local formation of a (screened) high-spin, and interorbital 
charge correlations drop abruptly, indicating the orbital 
decoupling [ ]. 

Some of these results reproduce previous ones 
[19, 35] and more generally the tendency to or- 
bital differentiation is indeed confirmed by a large 
amount of realistic simulations performed with sev- 
eral methods for treating the electronic correlations 
on top of DFT bandstructures: DMFT [40-47], varia- 
tional Montecarlo[48], Slave-spins [35, 49], Hartree-Fock 
mean-field [■">()], fluctuation-exchange approximation[51], 
Gutzwiller approximation [1 9]. 

All the studied stochiometric Iron-superconductors 
(representative of the '122', '1111' and '11' families) ap- 
pear to be in the orbitally decoupled bad metal state. In 
order to confirm that the scenario we have drawn in the 
main manuscript is general to iron superconductors, we 
performed calculations as a function of doping for each 




5.2 5.4 5.6 

(hole doping) 



5.8 6 6.2 

total filling (electron doping) 





1.5 




1.4 






s 


1.3 










Ig 


1.2 


c 





x 2 ^2 




- xz/yz — 













5.2 5.4 5.6 5.8 6 6.2 

(hole doping) total filling (electron doping) 

LaFeAsO - Fillings from 6.3 to 5.0 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 










x 




xz/yz 




xy - 




total filling = 6.0 • • • • 



1.1 1.2 1.3 

individual orbital filling 



1.4 



1.5 



Figure S5. Quasiparticle weights Z a and orbital populations 
n a (upper panels) as a function of doping for LaFeAsO. Z a as 
a function n a for each orbital a (lower panel). The striking 
linear behaviour signals the orbital-decoupling. 



Finally, it is worth pointing out that another, different 
theory for the differentiation of the correlation degree 
in iron pnictides, and based on the effect of spin-orbit 
interaction (that we neglect in this work) was proposed 
inRef. [521. 
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n a (upper panels) as a function of doping for FeSe. Z a as 
a function n a for each orbital a (lower panel). The striking 
linear behaviour signals the orbital-decoupling. 
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